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Abstract 

An integral representation of the partition function for general n- 
dimensional Ising models with nearest or non-nearest neighbours inter- 
actions is given. The representation is used to derive some properties of 
the partition function. An exact solution is given in a particular case. 
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As a simple prototype of a statistical mechanical system that undergoes a phase tran- 



studied in various ways. Since the exact solution of the free energy F and the sponta- 
neous magnetization for the two dimensional zero-field Ising model (on a square lattice) 
were obtained more than fifty years ago ||, many efforts have been made towards a 
detailed study of properties and the possible finding of exact solutions for the higher di- 
mensional Ising model or for a two-dimensional Ising model with nonzero magnetic field, 
for reviews see e.g., [4-7]. In this paper we study the Ising models by an "integration 
approach". We present an integral representation of the partition function for general 
n-dimensional Ising models with either nearest or non-nearest neighbours interactions. 
From this representation we prove that the partition function of n-dimensional homoge- 
neous (i.e. translation invariant) Ising systems on a square lattice is uniquely given by the 
eigenvalues of the related interaction coupling matrix. For one- dimensional homogeneous 
ferromagnetic systems with positive coupling coefficients the partition function satisfies a 
special equality which means that the variation of the interaction couplings of the system 
is related to the variations of the external magnetic field and the temperature. For some 
special cases of interaction couplings, in all dimensions, we get an exact representation 
of the partition function in terms of a Bessel function. We also give an approximate 
representation of the general partition function for the homogeneous Ising model, which 
shows that its leading term is given by the largest positive eigenvalue (resp. the largest 
absolute value of the negative eigenvalues) of the related interaction coupling matrix in 
the ferromagnetic (resp. antiferromagnetic) case. In the case of nearest neighbours in- 
teractions in two dimensions our integral representation is connected with the study of 
the Ising limit of the (j>\ double- well models, see e.g. ||. 

The partition function of the Ising models with the usual nearest neighbours interac- 
tions on a (^-dimensional lattice with A lattice sites is given by, see e.g. HJ, 



where <jj, i — 1,2, A, takes values +1 and —1, < ij > denotes the nearest-neighbour 
pairs on the lattice and K = J / kT, K' = H / kT with T the temperature, k the Boltzmann 
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sition for spatial dimensionality n > 1, the Ising(-Lenz) model M has been extensively 




(1) 



constant, H the external magnetic field and J the coupling constant. J is positive for 
ferromagnetic systems and negative for antiferromagnetic systems. 

We consider the partition function (|I|) in a more general form including both nearest 
and non-nearest neighbours interactions, 

Z A =Y J exp\KY J C%a t a 3 + K' £ a\ , (2) 

{<x s } I M=l i=l J 

where C° = (Cy) is a symmetric A x A matrix, 

C°=0, C°. = C°, i^j, z,j = l,2,...,A. (3) 

C?- stands for the interaction coupling between <7j and Oy Clearly if C° ij>0 = 1 and 
C?- = for all non-nearest neighbours pairs, the partition function @ is reduced to (]l|). 
In the following we call C° the interaction coupling matrix. 

Let < ij >i, < ij > 2 and < ij > a denote the next-to-nearest neighbours, next to 
next-to-nearest neighbours and a-th next-to-nearest neighbours interactions, respectively. 
In the following we consider systems such that all interaction coupling constants for any 
given a are the same, i.e., C <i j >a = c a , Wi,j = 1,2,..., A, a = 0,1, 2,..., A — 1. We also 
assume that the number of sites interacting with an arbitrary given site i is the same for 
all i = 1, A. We call these systems homogeneous. In particular this implies that we 
impose the periodic boundary condition on the boundary of the lattice if the lattice is 
a bounded portion of Z d (or equivalently the lattice is identified with a ci-dimensional 
torus) . 

Let li be the number of lines connected to a lattice site i. /j is then a topological 
invariant of the lattice and is called the topological link number associated with the 
lattice site i. For homogeneous systems with only nearest-neighbour interactions, the 
link number Zj is exactly equal to the number of the sites interacting with the lattice site 
i, and by homogeneity we have that Z 4 is independent of i, i.e., = I for all i — 1, A. 
For instance, for a periodic one-dimensional chain with nearest-neighbour interactions we 
have I = 2. For a two (resp. three) dimensional square lattice with nearest-neighbour 
interactions and periodic boundary conditions we have I = 4 (resp. 1 — 6). Therefore 
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for a homogeneous system in d dimensions with only nearest-neighbour interactions (i.e., 
such that C<jj> = 1 and Cg- = for all non-nearest neighbours pairs), the matrix C° 
has the following properties: 

ECJ = = Vj = l,2,...,A; (4) 

i=i i=i 

E Cg. = A/; Tr(C°) = 0. (5) 

In fact, a higher dimensional system with nearest-neighbour interactions can always 
be viewed as a lower dimensional system with non-nearest-neighbour interactions, and 
vice versa. For example, for a one dimensional chain with nearest and next-to-nearest 
neighbours interactions, 

1 2 3 4 5 6 
• • • • • • 

we have C° i+1 ^ 0, C° i+2 ^ 0, i.e., a lattice site % has interactions with the lattice sites i±l 
and i ± 2. This system is equivalent to a two dimensional band with nearest-neighbour 
interactions, 



2 4 6 8 10 




1 3 5 7 9 



where the double lines stand for nearest-neighbour interactions in the one dimensional 
chain and the single lines stand for next-to-nearest neighbours interactions in the one 
dimensional chain. If different interaction couplings in different directions are assumed 
then we have C° i+1 ^ C° i+2 . 

To describe the interaction properties for both systems with nearest and non-nearest- 
neighbour interactions, we define a generalized link number L of the lattice by, 

^ECS = EC& Vj = l,2,...,A. (6) 

i=i i=i 
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When the systems have only nearest-neighbour interactions L is equal to the topological 
link number I of the lattice. Instead of (|5|) we have, for homogeneous systems with either 
nearest-neighbour or non-nearest-neighbour interactions (in arbitrary dimension), 

E Cg = AL; Tr(C°) = 0. (7) 

Any real symmetric matrix M (over the complex numbers <C) can be diagonalized by 
an orthogonal similarity transformation. The matrix used to diagonalize M has as its 
columns an orthonormal set of eigenvectors for M. The resulting diagonal matrix has as 
its diagonal elements the eigenvalues of M. Let Ui be an orthonormal basis in € A (as a 
A-dimensional Hilbert space, with scalar product ■) of the (column) eigenvectors of C°, 
with eigenvalues A°, i — 1, 2, A, 

C°Ui = X°Ui, Hi ■ Uj = 5ij, i, j = 1,2, A, (8) 

(with Ui the adjoint vector to Ui). Let A be the orthogonal matrix that diagonalizes C°. 
Then 

A = («i,u 2 , ...,« A ), (9) 

AC°A = diag{\l,\l,...,\l}, (10) 

where diag{\1, A°, A° } denotes the A x A matrix having A?, A°, A° on the diagonal 
and elsewhere, and A is the adjoint of A. Moreover we have 

AA = 1, detA = detA = 1. (11) 

From (H) we also have 

= (12) 
Summing over the index i on both sides of the equation (0) and using (|J) we get 

(L - A° fc ) ][> lfc = 0. 

i=l 

Therefore if A° 7^ L for some = 1, A, then: 

A 

£A fc = o. (13) 



i=i 



Let us now consider the possibility that A° = L for some k, i.e., L is an eigenvalue of 
C°. We study the eigenvector x =column(xi, x 2 , ...,£a) £ (C A to the eigenvalue L (if it 

A 

exists) of C°. We have E C«^j = i = 1, 2, A. From @ we see that 

3=1 



I 1 \ 



Ul 



(14) 



V 1 / 



is an eigenvector of C° with eigenvalue L and with unit length (up to a global phase 
factor). By the orthonormality of the u\ we have U\U\ = 0, i — 2, 3, A. Hence from 
(|i~4]) the matrix that diagonalizes C° has the properties 

a r o, k ^ i 

E^= ^ (is) 

i=i [ vA, k = 1. 

Generally the matrix C° has both positive and negative eigenvalues. Let A+ aa; be 
the largest positive eigenvalue and A~ ax be the largest absolute value of the negative 
eigenvalues. Set A^ = ±A^ aa , ± e for any small constant e > 0. We set 



C + = X+I - C l 



(16) 



and 



C~ = C°- X~I, (17) 

where I is the A x A identity matrix. C^ 1 are then positive definite matrices with the 
properties: 



C± = ±\ ± , C% = C%, i^j, i,j = l,2,..,A, Tr(C ± ) = ±X ± A (18) 



and 



£qfc = £q± = ±A± TA vj = i,2,..,a, 



»=i 



8=1 



E C5 = (±A ± T L)A. 



(19) 
(20) 



The eigenvectors of C° are also eigenvectors of C 1 * 1 , 



C^^Afwi, wit/i Af = ±A ± ^A^, i = 1,2,..., A. 



± -r- \0 



(21) 



We have (with A as in (^)): 



AC ± A = diag{X£, Xf, ...,A±} 



(22) 



Corresponding to ([12]) we have 



(23) 



Set K + = J + /kT, with J + = J (resp. = J /fcT, with J = — J ) for ferromag- 
netic, i.e., with J > (resp. antiferromagnetic, i.e., with J < 0) systems. From (|T6|). 

A 

(|i"7[) and erf = A, we see that the partition function @ can be rewritten as 
i=i 

Zt = exp {±X ± K ± A} J] ea:p J -K ± £ Cj^- + AT' £ a t I , (24) 

R} [ ij=l i=l J 

where Z\ (resp. Zj) represents the partition function for ferromagnetic (resp. antifer- 
romagnetic) systems. In formula ([24] ) all the interaction couplings could be included in 
the elements of the matrix C ± . The parameter is no longer independent and could be 
scaled to be 1. Zj^ is obviously independent of e by the definition of X ± , (|T6|) and (|i~7|). 



Formula (|24|) can be transformed into an integration by using the following lemma, 
see e.g. @: 

[Lemma]. Let / and 5* be real-valued measurable functions on IR d such that F(X) = 
J M d f(x)exp {XS(x)} dx exists and S E C 2 (lR d ). For < A — > oo, F(X) is equal to 



-F(A) = (y) f E \detS"(x°)\^ [/(a; ) + 0(±)] {AS(x )} , 
provided S^) has a finite number of non-degenerate maxima x° such that 



(25) 



detS (x) Lo = (iet 



S'(x)\ x o = VS(x)\ x o = 0; 
d 2 S(x) 



dxtdxj 



^0, i,j = l,2,...,d 



and max xeM dS(x) = S(x°). The sum in (|25|) goes over all such maxima a; . 
Now let 

d 

S( x ) = - - I) 2 , XiGR, x = ...x d ) e iR d . 

8=1 



(26) 



Then the critical points x° of S are given by 
dS{x) 



Here we have 



dxi 



detS" = det 



A(x i - l)xi\ x o = 0, 



l,...,d. 



'd 2 s{x) 

dxidxj 



det(-(l2x 2 -A)5 ij ). 



(27) 



(28) 



From (|3) and (H) we see that all critical points (maxima) a; £ M d satisfying detS" ^ 
are given by x° = (x®), i = 1, d, for all possible combinations of values x° = ±1. We 
have max xeIR dS(x) = S(x°) = 0. 

Let C denote the set of all critical points. By the Lemma we get asymptotically for 
A — > oo, 



J Rd f{x)exp{-\ jjx 2 - 1) 2 | dx = (y)*8- 



y€C A 



for any / such that the integration exists. 
In particular we have 



lim \i [ f(x)exp\-X J2(x 2 - l) 2 ) dx = (2vr)f 8^ £ f(y) 



(29) 



We shall now apply this result to the case d = A, observing that X^ee f(y) amounts to 
sum over y with yi = ±1 (independently, for each 2 = 1, A) and identifying then yi 
with the Oi occurring in ( p4|) we see that the partition function (p4[) can be written as 



exp|±A ± K ± A}(^)t8 

I- J /.TV 



JR 



exp 



2tt' 

E Cg W% + K'J2 Vi -\ J2(y 2 - l) 2 j> rfy | 



(30) 



A— >00; 



1=1 1=1 

with the notation -F(A)|a^oo = hniA^oo -^(A). 

We call (p0| ) the integral representation of the partition functions Z^. Now we shall 
study its properties. First however we make a remark: 

Remark 1. By combining the first two terms in the exponential of the integrand in (|30|) 
into a quadratic form and making a translation of the integration variable y, — > yi — 
in Q3g), with 

(31) 



a ± = 



8 



we get 



exp{±\ ± K ± A}(^)h^ [ expl-K* £ C%{ Vi + a ± )( 2/j + a±) 
+K ± (a ± ) 2 £ -\J2(y- - l) 2 } dyl^oc 

ij=l i=l J 



+K ± (a ± ) 2 2 Cg-ADaw-a*) 2 -!) 2 U/U^x, 

A 

Therefore it is clear that the external magnetic field H = K'kT = — 2a ± J ± £Cg con- 



i=i 



tributes a global translation a of the critical points of the function (Eq) together with a 
f A 1 

constant factor exp < if ± (a ± ) 2 £ C^j > to the zero-field partition functions Z^(H = 0): 



Zt(H = 0) = exp{±A ± K ± A}(A)| 8 t f expl-K* £ C tjVm 

2tt J ei [ ij=l 



A 



(32) 



-A£(y 2 -i) 2 \dy\ 



A— »oo- 



i=l 



Conversely, a global translation of the critical points in the zero-field partition function 
is equivalent with the introduction of an external magnetic field. 

For convenience in Theorem 1 below we shall denote by A the number of points in 
a one dimensional sublattice in an n-dimensional square lattice, so that A n is the total 
number of points of the given n-dimensional lattice. From (|2]) it is obvious that for given 
K, K' and lattice number A, the partition function is uniquely determined by the 
interaction coupling matrix C°. Further we can prove the following: 

[Theorem 1]. For given K, K', the partition function Z\n of n-dimensional homoge- 
neous Ising systems (described by (0)) on a square lattice with A" lattice sites is uniquely 
given by the eigenvalues of the corresponding interaction coupling matrix. 

[Proof]. We first consider the one- dimensional homogeneous case. In this case the 



9 



interaction coupling matrix C° has the form 



/ ai a 2 
a A ai 



(33) 



V a 2 a 3 ■ • • ai / 

where a a , a = 1,2, ...,A are real constants representing the (a — l)-th order nearest 
neighbours interactions. The diagonal element a\ is in fact zero, a% = 0. Since the 
system is homogeneous we have 

A 

A even 

(34) 



a a — oa-q+2 * 



(\ = 2, 3, — , A even 



a = 2, 3, 



A + l 



A odd 



This implies that C° is symmetric. 

The matrix ([33]) is generated by cyclically permutating the elements of the former 
row to the right. We call C° a right cyclic matrix. By using the A x A permutation 
matrix P, 



/ 1 
1 





1 



\ 


1 




(35) 



the right cyclic matrix C° can be expressed as 



C° = ai + a 2 P + a 3 P 2 + ■■■ + a A P A " 1 . 



(36) 



The k-th component of the a-th eigenvector of the A x A permutation matrix P acting 
in C A is simply given by 

2vu 



exp 



| x (a-l)(&-l)j, A;,« = 1,2,...,A, 



(37) 



see e.g., JT^, [TT|. As the eigenvalues of a real symmetric matrix are real, both the real part 
and the imaginary part of the eigenvectors fl37D are also eigenvectors for the matrix P. 
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We take the k-th component of an eigenvector u a , a = 1, 2, A, of P to be normalized 
by the factor 1 / \/A as follows, 



(U a )k = < 



1 /27T, 

cosl— (&-!)(«-!)) , l<a</3, 



1 : sin (-T-{k - l)(a - 1) ) , A > a > (3, 



(38) 



I v^A VA 

where k = 1,2, A, p = (A + 1) /2 for A odd and (3 = (A + 2)/2 for A even. From (f5§) 
we see that «i is just the eigenvector (|T4]). 

The orthogonal matrix that diagonalizes P is then given by 



A = (u h u 2 , ...,ma)- 



(39) 



It is direct to check that the matrix A obtained in this way satisfies the equations (O) 
and fll5|) . From (|36|) , (|i"7|) and fll~6f) , the matrix A also diagonalizes C 1 * 1 = ±X^I =F C°, in 
the way given by formula 



Changing the integration variable to be ^ Aj?/? in the integration (|30| ) and using 

i=i 



the relations (pTTf) and (p2|) we get, 



A A 

+if ^ - A 51 (AjAkyjVk - 1) 2 } 

ij=l ijk=l 



(40) 



From (|38|) and (|39|) , we see that A = (A^-) is independent of the elements of C^. Therefore 
for one dimensional homogeneous systems the partition function (|40D only depends the 
eigenvalues Xf , i.e., the eigenvalues of the interaction coupling matrix C°. 

For two dimensional homogeneous square lattice systems with A 2 lattice sites, the 
interaction coupling matrix C°, denoted here by C\, may have various forms according 
to the ways of numbering the lattice sites. We number the lattice sites from left to right 
and from the first line to the last line of the lattice. The matrix C® is then of the form, 

/ C° a 2 I •■■ a A I \ 
a^I C° ■ ■ ■ a\-\I 



C° 2 



C Q ®I + I®C Z 



(41) 



V a 2 I a 3 I ■■■ C° J 
11 



where C° and a a , a = 2, A, as in fl33|), and / is the A x A identity matrix. 

C 2 is again a right cyclic matrix generated by permuting the matrix blocks (C°,a 2 I, 
(I3I ,...,a&I). The eigenvectors of the interaction coupling matrix C 2 have the form, 

v a p = u a ®u p , a, (3 = 1,2, A, (42) 

with u a , up as in (|38|). The matrix that diagonalizes C 2 is then given by 

A = {v n ,v 12 , v 21 ,v 22 , «aa)- ( 43 ) 

A is then a A 2 x A 2 othonormal matrix which is independent of a 2 , ...,Oa- Obviously 
the matrix (f43|) also diagonalizes the positive definite matrix = ±X ± I ® I C 2 
in the integral representation of the partition functions. Therefore for two-dimensional 
homogeneous systems the partition functions are also only determined by the eigenvalues 
of the interaction coupling matrix for given K, K' and A. 

Similarly, for an n(> 2) -dimensional homogeneous system, the interaction coupling 
matrix C° has the form, 

/ C°_! a 2 I ■■■ a A I \ 

-1 ' " ' a A-l-< 



(44) 



V a 2 I a 3 I ■■■ Cl_ x ) 
where C°_ x is the interaction coupling matrix of an (n — l)-dimensional homogeneous 
system. 

The matrix that diagonalizes C° and = ±\ ± I n =fC°, with I n the A™ x A™ identity 
matrix, is given by 

A = (u ai ® u a2 (g ... <g> «„„), cKj = 1, 2, A, Vi = 1,2, ...,n, (45) 

where w ai is given by ([38]). This proves the theorem. ■ 

[Theorem 2]. For one- dimensional homogeneous ferromagnetic Ising systems with in- 
teraction coupling matrix C° given by fl3"5|) and coupling coefficients a, > 0, i = 2, 3, A, 
the partition function Z\ satisfies 



A xy+ 



i=2 



+ rpdZ^ + R dZi 



dT OH 



dZ 



4 =2 9a i 



K + A(A — \)LZ\. 



(46) 



12 



[Proof]. From the partition function of one-dimensional ferromagnetic homoge- 
neous systems takes the form, 



exp{\ + K + A}(±)h$ f^exp^-K+Y.Xfyl 

A A 

+K' £ AjVj - A £ (AjAkVjVk - l) 2 \ dy\\^oo 

ij=l ijk=l ) 



(47) 



B + / exp\v + )dy\ 

J]R A 1 > 



A— >oo> 



where A„-, i,j = 1,2,..., A are the elements of the matrix 



\+ _ \+ _ \0 _ \+ i_ 



e — A°, A° is an eigenvalue of C°. For simplicity we have set in Q4"7| ) 



5+ = esp {A+^+A} (— )*8 



2tt' 



(48) 



r+ = - K + J2 Kv" + Am - a E (A,A ife%m - 1) 2 . (49) 

ij=l ijk=l 



1=1 



By ( j24|) and the definitions ([TBI) , (|i~7D of C + , as we stated before, the partition function 
Z\ is independent of the positive real number e. Hence we have 

dZt 







(50) 



On the other hand using Lebsgue dominated convergence one shows that the application 
of d/de to the right hand side of fl4"T|) yields 

K+KZl + exp \X + K+A\ (— )*8*- 



. f A A A 

/ a exp -If + £ Xfyf + K'Y, Am ~ A £ (AjAkVjVk - I) 2 

J7R [ i=l tj=l ijfc=l J 



(51) 



1=1 



A— >oo> 



(In fact the derivatives of Z\ with respective to e and, in the following, to T, H and dj, 
can be exchanged with the integration and the limitation in (fl"7D because from ([24]) these 
derivatives only change the form of the function / in ( p9|) ). Equating ([H]) with (|50|) we 
get ^ 

AZ+-S+ / exp{r + }E2/ 2 ^|A^oo = 0. (52) 
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In a similar way we get from ([17] 
dZi 



OH 



kT 



ij=l 



dy\ 



A— >oo- 



(53) 



and 



d_Zl 
dT 



x+A (-£) z > 



J+ A H A 

1=1 ij = l 



(54) 



dy\ 



A— >oc 



From (|53|) and (p4|) we deduce that 



dZ+ + H dZt 



dT T OH 



Using (fj^) one gets 



dZt TT dZt 
T—^ + H- ' 



dT OH 



+ B+J MA exp{T+} 



B+ L exp i T+ } 



~T 



i=l 



(55) 



dy\ 



A— >oo- 



i=l 



A— »oo- 



(56) 



(57) 



Now we discuss the properties of the eigenvalues of the matrix C° given by (p3[). We 
first recall a general result of matrix theory. Let g(x) be polynomial in x of degree I with 
roots (pi,p 2 , -,/>/), 

#(a;) = cio^' + d\x l ~ l + ... + di^ix + d\ 

= (-l)'do(pi - x)(p2 - x)...{pi - x). 
Let Xi, i = 1,2, ...,m, be the eigenvalues of a given m x m matrix B over C, acting in 
C m . Then 

|p/- J B| = (p-A 1 )(p-A 2 )...(p-A m ). 
Replacing x by £> in (|57|) we get 

<?(B) = (-l) l d (pil - B){p 2 I - B)...( Pl I - B). 



Hence 



\g(B)\ = (-ir^\( Pl I-B)\\(p 2 I-B)\...\( Pl I-B)\ 



l m 



■i) mi cnn(ft- A .)=n^)- 



(58) 



k=l i=l 
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see e.g. ||TT] . 

From (|36| ) we see that C° is a polynomial in the A x A permutation matrix P acting 
in C A . As \pl — P\ = p A — 1, the eigenvalues of P are 

'2m(a- 1)1 



A 



a = 1,2, A. 



(59) 



Noting that all diagonal elements of C° are equal to a±, from (|58|) we have 
|C° — A°7| = \a l -X° + a 2 P + a 3 P 2 + ... + a A P A - 1 \ 

= f[[ai-X° + a 2 ( a + a 3 ( 2a + ... + a A C (A - 1)a 

a=l 

The eigenvalues of C° are then given by 

A° = a 2 ( a + a 3 ( 2a + ... + a A C (A - 1)Q , a = 1, 2, .., A. 



(60) 



(61) 



By the definition (|6|) we see that A? = Sa=2 a « * s °f special significance. It is just the 
generalized link number of the lattice, L = A°- For our case, a a > 0, a = 2, 3, A, and 
the link number L is also the largest eigenvalue of the interaction coupling matrix C°, 



a; 



A? = L, 



which means Xf = A+ aa . + e — A? = e and hence 

dXt 



<9A° 



0. 



1,2,..., A. 



From ( |6"T| ) A° satisfies 



A <9A° 

7=2 



7 



A - 1, a = 1, 
-1, 

The eigenvector corresponding to A° is w a given by (p8j) . 



u « - 2^ Q 



C% = X a u a . 



(62) 

(63) 

(64) 
(65) 

(66) 



From fl47|), (@) and (H) we have 



E- 

8=2 



9a,- 



n A 

EE a < 



9z+ <9a° _ A ^ 



i=2 M =l ^ M 



a<9A° 



LK + KZl + 5+ 



^ + E x °iV 2 



i=2 



(67) 



A— >oo; 
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where the relation (p2]) has been used. 
By ((47|), (§§ and (§3|) we also have 



E 

j=2 



EE 



E 



' <9A° ~ M 

At=l M 



i^+A(A - 1)^ -B + I exp {r+} 



^ + EW 



i=2 



(68) 



A— >oo- 



From (|56T), (167T) one has 



i=l 



da; 



dT 



dH 



LK+AZi 



B+ (ex P {T + }(-K + Lyl)dy\ 



(69) 



A^oo- 



By fl52|) and (0) one gets 
A dZt 



E 

i=2 



K + A(A - 2)Zt + K + B+ f exp {T + \ y\dy\ 

JM A 1 > 



A— >oo- 



(70) 



Equations (|69"D and (|70|) then give 



A dZt dZ 

2^ a i— ^ T 



+ 



9a. 



<9T 



<9# 



9a, 



1=2 ^"t ^ j=2 

= LK+AZl + LK + A(A - 2)Z^ = LK+AZ^(A - 1), 



which is just the equality ([46]). ■ 

Remark 2. Equality ( f4~6|) means that the variation of the interaction couplings of a 
system is related to the variations of the external magnetic field and the temperature. 
As the free energy F + , internal energy u + and magnetization M + per site are given 
respectively by 



u 



-T 

M + 



F + = --kT log Z + , 
A 

^ (.± ]ogZ A- T2kd 



dT \ A J A dT 

8F+ Tk d 



\ogZ- 



3) can be rewritten as 



dH A dH 



3F+ 



\ogZ~ 



u 



+ M + H - J> 4 + L)-^— = J+(A - 1)L. 



i=2 



(71) 
(72) 

(73) 
(74) 
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([74]) gives the relation among the internal energy, magnetization, external field, free 
energy and interaction couplings of one-dimensional homogeneous Ising ferromagnetic 
systems. 

[Theorem 3]. For one-dimensional homogeneous Ising ferromagnetic systems with in- 
teraction coupling matrix C° given by (|3~3"D and coupling coefficients a« > 0, i = 2, 3, A, 
the correlation function between spins i and j satisfies 

1=2 1 1=2 1 

[Proof]. Let f(a) be an entire function of c^, i = 1,2, A, independent of a«, T and 
H. The average < f(a) > of f(a) is by definition given by 

< f(a) >= ±- £ f(v)exp I K £ Cg^ + AT'X> i l, (76) 

^ A {o-i} I »J=1 »=1 J 

where Za is given by (fj). 

For ferromagnetic systems, using our integral representation we see that we can write 
(|76|) in the form: 

< m >= J rA f(a)exp {r + } dy\ x ^ = (77) 

where Z~j[, B + and Y + are given by (]47|), (|48| ) and ( f49| ) respectively, and 

Z+(/) = E+ £ a /(<r)ezp {r + } dy]^. (78) 

Let (7 be a linear operator defined on smooth functions of a{, T and H by 

A 8 d 3 A 8 

v = X«eZ + T Wr + H M + L X*;- (79) 

From Theorem 2 we see that 

UZt = K + A(A - \)LZ\. (80) 
It is also straightforward to check that 

UZ+(f) = K + A(A-l)LZ+(f). (81) 
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By ©> © and (0) we get 

</(*)> = ^uzt(f) - ?UQuzt 

i z+(f) (82) 

= ^ + A(A-l)LZ+(/)--A^i^+A(A-l)i:Z+ = 0. 

Z A \ Z 1) 

For the correlation function =< (TjOj > — < (Xj >< <7j >, we get using the result (|82|): 

Ugij = U < OiOj > —U(< (Ti >) < (Jj> — < CTi > U(< (Tj >) = 0, 

which is just the formula (|75|). ■ 

By using our integral approach to the Ising model, we have discussed some properties 
of the partition function for a one-dimensional homogeneous case. We now remark that 
for some zero magnetic field cases, in arbitrary dimensions, the partition function can be 
given exactly in terms of special functions. 

Let A be the matrix that diagonalizes C as in (p2|) . Let us consider the Ising system 
described by (f24|). 

[Theorem 4]. If is such that A satisfies 

A 

^ ] AijAi^Ai m A{ n 8 jk^km^mn^n (83) 
i=l 

for some constants b n , Reb n > 0, n = 1, 2, A, then 



n 



K(Xf-2X) Ki (K\Xf-2Xr 



(84) 



A— > + CXD? 



l = \\\ Xh ±\ 8Xb t 

where K\ is the Bessel function of imaginary argument equal to \. 

A 

[Proof]. The partition function of zero-field is given by (^). Set y,i = AjjPj in 
(|32|), A = (Aij) as in (^). By using (^), ([TT|) and fl83|) we have 

E C±y iVj = Y2Xtpl (85) 

AAA A A A 

Y,Vi = E(EA,^) 2 = E J2 A a A mPk = E? 2 ( 86 ) 

j=l j=l j'=l J,fc=l i=l i=l 
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and 

AAA A 

J2yt = J2 A ij A ik A im A in p j p k p m p n = J2biPt- (87) 

7=1 7=1 j, k, 771, 71=1 7=1 

A 

Substituting = /^A^pj, ([S5l), (|S6|) and (|87D into (^) we get 

3=1 



Z£(H = 0) = exp{(±X ± K ± -A)A}(A)M 

Z7T 



V exp (-#± £(Af - 2AK 2 - A £ 6^1 dp| 

I 7 = 1 8=1 J 



A^oo- 



From the formula (see e.g. fll 



L exp i-^ - 2vx2 ) dx = ivf exp {h} K * ih) ' Refi > °' 



(88) 



we obtain (R^). 



In the following we present an asymptotic formula for the partition function (0) in 
the general case. We first make a Laplace transformation to the integral representation 
of the partition function ([30]) . 

From the formula 



for any symmetric strictly positive definite matrix B, (|30|) can be reexpressed as 



Zf = c±A2 



J1R JR { i,j=l M=l 



A 



A— >oo> 



i=l 

where a 1 * 1 as in (HTft and 



c ± = exp 



|±A ± K ± a| (— )Mexp I 7T ± (a ± ) 2 V C& \ (—)^VdrtC±. (90) 
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By a translation yi — > yi + 1 and a rescaling ?/j — > yij y/X we get 



Z± = c±A^ 



L L e ^ I - K± s + 2iK± e ( a± + 

JR JM { i,j=l i,j=l 



i,j=l i=l 



(91) 



/ A / a ex p I ~ K± E + 2 ^ ± E ( a± + 

E^% - E(y + ^ + 4y?)| dpdyU-.cc. 

[Theorem 5]. For given i^, if', up to order 0(1/A 2 ), the partition function of any 
n-dimensional homogeneous Ising system is given in terms of the link number L and 
the largest positive eigenvalue A + (resp. the largest absolute value A~ of the negative 
eigenvalues) of the related interaction coupling matrix for ferromagnetic (resp. antifer- 
romagnetic) systems, according to the formula 

Z£ = exp{{±K ± L + K')A}\ [ l + j[^ + ^{2K ± {±\ ± T L)-K') 

+ x i ( 92 ) 



+ -{2K±{±\± T L)-K'Yt—\ ± 



[Proof]. By expanding the integrand in ( |9lD to order — , before taking the limit 

A 

A — > oo, and using the integration formulae 

J^x 2n exp{-px 2 }dx = ^~ p) l )U ^ P>0, (93) 
where (2n - 1)!! = 1 • 3 • 5 • ... • (2n - 1), and 

/CO 
x 2n+1 exp{-px 2 }dx = 0, p > 0, (94) 
-co 
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we have 



JM A JM A 



exp{-4yf)} 
1 



exp | -K± J2 C %ViVi + 2iK± J2 (« ± + WtjPi \ 
-Lf-4(jf±)» £ CgcSftp^M + 16 £ 

\ i,j,k,l=l i,i=l 



dpdy + 0(—) 



= c± [ [ expl-K+Yl C %ViV 3 + 2^ E (a* + l)C±p, 1 exp {-4^ 2 

JK JR { i,j=l i,j=l ) 

+16$>?-16iir ± E 

*=i ij=i y 



dpdy + 0(^ 



= c± 



,/n\ A f 1 / 3 . 3^ A ^ 



2 y [ A ^16 8 ^ 

(^±)2 A 



ijPi 



E cgcgpipj 



dp + 0(-) 



i,j,k=l ) J 

(^Pj / kA j-jf* EcgCft - *( a± + i))(Pi - *( G± + 
^ ± (« ± + i) 2 E^ 



i)) 



1 / 3 . 3iK± , 
1 + A I6 A "— 



ij',fc=l / 

where it has been taken into account that 



dp + 0(—), 



O p ,y{-^) standing for the remainder in the expansion of the integrand in (plf) up to 



1 

A 5 ' 
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Let pi — > pi + i(a ± + 1). 



7 ± 



;± (^) Ae ^{-^ ±(o±+1)2 ^ 1 ^} L exp { 



A 



i 



E 05^ + ^ + 1))^ + 1(^ + 1))) 

i,j,k=l J 

K±{*± + l) 2 £ Cjl / A exp ( E CyftPi 

ij=l J JK I ij=l 



1 + i ^^.Ec^^^^i^ 



i,j,k=l 



dp + 0(-), 



where 



JJR I ij=l J ij,*=l 

= / a exp i ~ J ^ ± E A i t ^ 2 f E ((C ,± ) 2 )ii^m^jng m gndg 

JK I i=l J i,j,m,n=l 

= [ A exp \-K± E E^) 2 ^ 2 ^ 
JK I i=i J i=i 

- 1 A~ K± t AW ) s <A - )2? ^ 

Therefore we get 



TrC = 



exp |±A ± iT ± A - i^ ± (2a ± + 1) E j 

ij,k=i 



1 + M^A+^E^ 



A I 16 



+ o( F ). 



8 
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Using (|19|) and (pO[) we obtain 



Z 



exp{{±K ± L + K')A} 1 + 



+ — (2K ± (±\ ± tL) -K') 2 t 



A / 3 
A 1 16 " 



+ —(2K ± (±X ± tL)-K') 



(96) 



8 



From (pBD one gets a corresponding asymptotic representation of free energy of the 



system, 



F 



-kT- log Z A = -A;T(±K ± L + if') 




(97) 



Conclusions From the integral representation of the partition function for general 
n-dimensional Ising models with both nearest and non-nearest neighbours interactions we 
have proved that for given K, K', the partition function of n-dimensional homogeneous 
systems on a square lattice is uniquely given by the eigenvalues of the related interaction 
coupling matrix. For one- dimensional homogeneous ferromagnetic systems with positive 
coupling coefficients, the partition function satisfies a special equality which means that 
the variation of the interaction couplings of a system is related to the variations of the 
external magnetic field and the temperature. For some special cases of interaction cou- 
pling, we obtained, for a Ising model in n-dimensions, an exact solution for the partition 
function in terms of Bessel functions. We also calculated the partition function of the 
n-dimensional Ising model to order 0(1/ X 2 ). It turned out that the leading terms of the 
partition function for any homogeneous systems in arbitrary dimensions are given by the 
largest positive eigenvalue (resp. the largest absolute value of the negative eigenvalues) 
of the related interaction coupling matrix for ferromagnetic (resp. antiferromagnetic) 
systems. 

The advantage of our integration approach is that one can analyse the partition func- 
tion of the Ising model for various interaction couplings on lattices of arbitrary dimensions 
in terms of integrals. More results could be obtained by studying the properties of other 
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interaction coupling matrices. It would also be interesting to study the relation between 
the largest eigenvalues of the interaction coupling matrices in our integration approach 
and the largest eigenvalues of the usual transfer matrix approach to Ising models, see e.g., 
0, ^ . Moreover the integral representation formula can also serve as an alternative way 
to study the correlation functions, and as a mean to provide numerical approximations 
to the thermodynamic functions of the Ising model. 

ACKNOWLEDGEMENTS: The original idea for the integral representation was 
formed in discussion of the first author with the late Raphael H0egh- Krohn, in the 
fall of '87. We dedicate this work to his beloved memory. 
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